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Abstract 

CN ' We study the t-J model on a ladder using the slave-fermion-CP 1 formalism which 

o : 

has been used successfully in studying lightly doped high-Tc cuprates. Special at- 



tention is paid to the dynamics of composite gauge fields and the natures of quasi- 
excitations. The slave-fermion-CP 1 approach explains many aspects of the ladder 
system observed by experiments and numerical studies in a natural and coherent 



manner. We first obtain the low-energy effective model by integrating out half of 

o 
o 



the CP variables (the Schwinger bosons) assuming a short-range antiferromagnetic 



order (SRAFO). The spin part of the effective model is the relativistic CP 1 model. In 
^ . the single-chain case, there appears a topological #-term with 9 = tt, as is well known. 

H : 

On the other hand, in the two- leg ladder case, we have 6 = 2tt. The dynamics of the 
composite gauge boson strongly depends on the value of the coefficient of the #-term. 
This fact explains why the quasi-excitations in a chain and ladders of even legs are 
different. For a ladder, the gauge dynamics realizes in the confinement phase, so the 
quasi-exciations are charge-neutral objects like spin triplet and electrons, etc. The 
effective model reveals attractive force between holes, which generates superconduc- 
tivity in SRAFO. The symmetry of hole-pair condensation should be of the d-wave 
type. 
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1 Introduction 



In the last few years, quasi-one- dimensional systems of strongly correlated electrons 
are of great theoretical and experimental interest in condensed matter physics]!]]. 
Studies on these systems are expected to shed light on the mechanism of high-Tc su- 
perconductivity. At present, it is known that spin excitations in an antiferromagnetic 
(AF) Heisenberg model put on an even-leg ladder have an energy gap, though those 
on an odd-leg ladder are gapless. Very recently, superconductivity has been observed 
in (weakly interacting) doped spin ladders with a spin gap@, and the symmetry of the 
order parameter seems to be of a d-wave type as some theoretical studies predict @. 

The metallic phase of high-T^ cuprates is anormalous, and a phenomenon called 
charge-spin separation (CSS) is expected to occur J|]. In mean-field theories (MFT) 
of the t-J model in slave-particle representations, the CSS is implemented naturally. 
However, the phase degrees of freedom of MF's behave like (composite) gauge fields 0, 
so careful study is necessary in order to check the reliability of MFT. In the previous 
papers H [7j, we argued that the CSS takes place at sufficiently low temperatures 
(T) in the two-dimensional (2D) t-J model. To derive the result, a gauge-theoretical 
treatment was essential. 

In the present paper, we shall study the t-J model on a ladder of two legs by 
using the slave-fermion (SF) representation. The gauge-theoretical method is also 
applicable for this model straightfowardly. It clarifies why the quasi-excitations in 
chain and ladder systems are different so much; In a chain, the full CSS takes place 
and holons and spinons are quasi-excitations, while in a two-leg ladder, we shall see 
that the system is in a confinement phase and the gauge-charge-neutral objects like 
spin triplet and bound states of holons and spinons are low-energy excitations. This 
result explains the numerical studies in Ref. ||. 

There are field-theoretical studies for the ladder systems, most of which use the 
bosonization techniques ||. We want to stress here that the gauge-theoretical study 
is useful not only for the 2D systems of strongly-correlated electrons but also for the 
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quasi-one-dimensional systems. Furthermore, we believe that it gives us an universal 
and coherent understanding of a wide variety of strongly-correlated electron systems 
including the fractional quantum Hall effect. 

This paper is organized as follows. In Sect. 2, we shall introduce the t-J model 
in the SF-CP 1 representation. This procedure is closely related with our previous 
work which studied the 2D t-J model|jlO[. In Sect. 3, we shall obtain a low-energy 



effective model by integrating out half of the spin variables (CP 1 variables) assuming 
a short-range AF order (SRAFO). The effective model, as a result of the existence 
of SRAFO, contains important interactions among smooth spin variables and hole 
field. In Sect. 4, we shall study the spin dynamics in the chain and ladder t-J model. 
The continuum limit of the spin part of the effective model is the CP 1 model with a 
topological term. Coefficient of the topological term depends on the number of legs 
and strongly influences the dynamics of composite gauge bosons. As we explained 
previously f|, |?J, the mechanism of CSS in strongly-correlated electron systems is 
described by a (de) confinement phenomenon of composite gauge fields. In Sect. 5 we 
shall study the dynamics of composite gauge bosons on a two-leg ladder using the 
method presented in M and identify quasi-excitations. We find that the system is in 
the confinement phase down to T = 0. So the quasi-exciations should be gauge-charge- 
neutral bound states of holons and spinons, like a spin triplet (magnons) and electrons. 
In Sect. 6, we shall study the possibility of superconductivity and its symmetry. The 
effective model involves an attractive force between holes in the SRAFO background 
which enhances a d-wave hole-pair condensation. In Sect. 7 we present discussion. 

2 The t-J model in the SF-CP 1 representation 

We consider the t-J model on a two-leg ladder whose Hamiltonian is given by 
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-» -» 1 -► -» 1 

+ H [ J H (Si+l,aS i>a ~ jn i+ i ia n ita ) + J'(S itl S it 2 - jTli^Tlifl} . (2.1) 

t 0=1 4 4 

Here the suffix a(= 1, 2) distinguishes one of two legs, i labels sites along each leg. 

is the electron operator with the spin a(= 1,2), and Si ja and n.i jQ are the spin 
and number operators on the site (i, a). The physical states must satisfy 

rii, a < 2, (2.2) 

on each site. In the SF formalism, the electron operator is expressed in terms of the 
bosonic spinon operators cii^ a and the fermionic holon %p iya operators as 

Ci,a,o- = ll>i j a a i,a,cn (2-3) 

and the physical state condition (|2.2| ) becomes 

Yl a i,a,* a i,a,a + ^{a^i,a = 1- (2-4) 

As in the previous paper fllP ], we solve the condition ( p.4|) by rewriting the operator 



a^a in terms of CP 1 operator z^ a as 

a- = (1 - t/jT ?/)• V^z- 

= (1 - ^la&i^Ziw (2.5) 

5Z4,o,«7*W = 1- (2-6) 

This representation is quite useful especially for the lightly doped case of the t- J model 
and the AF Heisenberg model 0, |TT| . We shall treat and z^ a ,<j as fundamental 



variables and employ the path-integral formalism (See Ref.[[T(J for detailed discussions 
on this formalism). 

From (|2.3|) and (|2.5|) , it is obvious that there appears a local gauge symmetry in 
the SF-CP 1 t-J model; 

^ e^.«-0i,«. (2.7) 



The electron operator Ci ta , a is invariant under this transformation since it is a com- 
posite of a spinon and an anti-holon. One may expect that there appears composite 
gauge bosons with respect to fl2.7Q . As we showed in the previous papers, these 
gauge bosons are introduced as auxiliary fields in path-integrals, and their dynamics 
specifies the nature of quasi-excitations. 

3 Short-range AF and low-energy effective model 

In this section we shall obtain the low-energy effective model by integrationg out 
half of the CP 1 variables Zi^ a (e.g., the CP 1 variables on the odd sites) assuming a 
SRAFO. 

In the path-integral formalism, the partition function Z is given as 
Z = J [Dz\[D4)\ exp [ J drA(r)] , 

A ( T ) = -^(^Zi,a,aZi,a,a + $i,a>Pi,a) - H. (3.1) 



The J-terms in H is explicitly given in terms of Zi^ a and ipi ja as 

i,a 

J' 

Aj> = ~irJ2Pi,lPiA( Z i,l Z h2)( z i,2Zi,l) ~ 1], (3.2) 
A i 

where p\ a = 1 - i>i >a il>i,a- 

From ( 13.21) , it is obvious that the SRAFO configurations, Sj ~ —Si for a nearest- 
neighbor (NN) pair give dominant contributions to the path integral. One can 
express a spin Si and its time- reversed spin — Si as 

— * 

Si — Zi<TZi 

Si — Zi(JZi 

z a = taa'Za', ^12 = — ^21 = 1, (3-3) 



5 



where a is the Pauli spin matrices. Therefore the assumption of SRAFO is written 

as 

Zi ~ Zj, = NN sites. (3.4) 

To integrate out half of the CP 1 variables, say, those on the odd sites, around the 
SRAF configurations (|3.4|) , we pair every odd site with a NN even site and denote 
the CP 1 variables in a pair as z Q and z e . For J > J', we take z e = Zi + \ A for odd-site 
z = Zt,ai whereas for J' > J, z e = z^a for z Q = z i>a where a = 2 for a = 1 and a — 1 
for a = 2 (see Fig.l). 

The odd site CP 1 variable z Q is parameterized as follows by using its reference 
coordinate z e , 

Zo = PoeZe + QoeZe, Qoe = (1 ~ PoePoe)^ U oe , U oe G U(l). (3.5) 

We substitute ( |3.5| ) into A(r) and expand it in powers of p oe up to 0{p 2 oe ); 

A = A + A p + O(p 3 oe ), (3.6) 



A = Yl [ - V'eV'e - VW>o + (Po - P\)ZZ + fi c {pl + p\) 
o£odd 

+ \p 2 o J uP 2 ou {.ZZ ov ){z ou z) - Y t uPoPoubPoi ! ouU oe {ZouZ) 

v=u,d,s u=u,d,s 

+Uoe{zz ov )i) ov i) ] , (3.7) 

where the suffix u(u, d, s) denotes the three NN directions from an odd site o, and 
we omit the suffix e as z = z e . fi c is the chemical potential to enforce the hole 
concentration to be (ipiipi) = 8. J u ,t u denote 



, J, z/=u,d 
J„={ (3.8) 



J', u=s 



t, z/=u,d 

(3.9) 

t' z/=s. 
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A p is given by 



Ap — 2^ ( PoeM a p oe + p oe k oe + l oe p oe )i 
o£odd 



1 — >• 1<— 

2^o^r - 2 d rpl + 2 pl ZZ 

+ \ y }l J vPlp 2 oA(ZovZ){zZ ou ) - (z ou z)(zZ ov )} 



1 



-p 2 U oe zz - ^ \-^- P IpI v {zz ov ){z ov z)U 



+tvPoPou1pov*Po(zZ, 

2 i ■ - ■ r 



toe Po^oeZZ 

+t v p pou'4>oi> ov [Zqi/Z ) . 
We do Gaussian integration over p oe 's; 



ododd 



At low T,T <2J+ J' this reduces to 10 



/drAi(r) = / ' dTdT r £loe(T)M- 1 (T,T')k oe (T') 
J J 

2 r 

~ 2 J + J' y rfr^/ oe (r)A; oe (r). 



(3.10) 



[ll>oi>auUae(z 01/ z) + U oe (zZ ov )lp ov %l) } 1 (3.11) 



(3.12) 



/" [dp]exp [ /" rfrAp] = Yl (detM ) _1 • exp / drAiir) . (3.13) 



(3.14) 



From ( |3.7| ) and (|3.14j) , the relevant terms of z and t/> are readily obtained. For the 
CP 1 field z, 



A* 



A\ = 2^(2J + J')- X [(zz)(zz) 

o 

J V J'- 

+ ^Z-^-^{ zz ou){z ov z){zZ ouI ){zq v iz) 
vv 1 

+( zz ) X -£(z z ov)(ZovZ) + (ZZ) X -^{ZZ 0V ){Z 0V Z) 



(3.15) 



Similarly, for the hole field tp, we obtain 



K 



K + T + Ti + T 2 , 

— ^2 \$(Dt + m)ip + fj(D T — m)r] 

o 

mi JufjoVo^ov1pou{zZ v)(z ol/ z), 
a v 

tuiK'tpouVo + c v ri i\) ov ), 

O V 

{zz ou ) 



m = jj c + 2 J + J', 



T 2 



-{zz oy ) + 
-{z ov z) + 



2J + J' 

(z ou z) 
2 J + J' 



2(zz) + Y Jv>{z ov iz){zz ov i) 

v' 

2(zz) + Y Jv'{zz ou ,)(z ou ,z) 



2 J + J 



7 tutu'11o^ov^ou'Vo{ZouZ){zZ ou i), 



o vv 1 



where 



(3.16) 



(3.17) 



D T = d T + iA T = d T — (zz). 
Above, we have defined 

Vo = U oe tp , (3.18) 

which transforms r] — > e l ^ e T] under ( 2/7) . The action of the effective lattice model is 
thus given as A eS = A z + A\ + Af el . 

T and T 2 in the effective model show that there appear effective interactions 
between holes as a result of the SRAFO. One can expect that a superconducting 
phase appears when weak-inter-ladder interactions are included as recently observed 
by experiments 0. The order parameter for superconductivity is the following hole- 
pair field, 

M ou = 4> ou (zz ou )r] , (3.19) 
which has the electric charge +2e and invariant under (12. 
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4 Spin dynamics in chain and ladder 

In this section, we study the spin part of A e s by taking the continuum limit. An 
essential difference appears between the ladder and chain. 

We shall first consider the case J' < J and focus on the imaginary term in the 
effective model, i.e., the last two terms in ( 3 .151 ). 



I z = Tpj—J, J » [( zz ) {zz ov ){z ou z) + (zz) {zz ov )(z ou z) 

= 2 J + J' ^ J -[C^)(z ou z) + (tz)(zz ov )\ + ... (4.1) 

We introduce (continuous) coordinate x = ai (a is the lattice spacing). For smooth 
configurations of z, we obtain 

JdrF ~ —L— a ^2(2J + J') jdr(p^D T z-D^D x z) 
~ J dxdr{D x zD T z — D T zD x z^j 

= 2mQ, (4.2) 

where = + iA^ = — (zd^z), and Q is the topological charge which takes 
integer values for smooth configurations of z since it is a wrapping number from S2 
of the x — t space to 0(3) of the S space. On the other hand, for a single quantum 
spin chain, 

r dTP cham = <KiQ. (4.3) 



The reason why the coefficient of the topological term doubles in the ladder case is 
simply because there are twice as many degrees of freedom per unit length in the 
ladder case. 

Before going into the detailed study on the effect of the topological term, let 
us obtain the continuum limit of the remaining terms in and A\ in (3.15). By 
straightforward calculation, we obtain 

^ o 



~ -(2 J + J')a J ' dxD x zD x z, 
A\ ~ — - f dx\- a- v D^z~D T z + {J + J' /V2) 2 a~D x ~z~D x z\ + F . (4.4) 

Therefore, the continuum limit of the spin part of A z = Aq + A\ is the relativistic 
CP 1 model with the topological ^-term. 



A CP = — [ dx V D^zD^z + I 



(J,=T,X 

., 1 + J'/2J 



gf = 1 ' =, (4-5) 

A + (2 - y/2)J'/J 



where we have rescaled the imaginary time as r — > v z r and t> z = [J( J+(2— \/2) J')] 1 / 2 a 
is the "speed of light" of the present system. Recently, the 0(3) nonlinear-a model 
with 0(3) variables n(x) (n(x) -n(x) = 1) was derived as an effective-low energy model 
of the AF Heiseberg models on a ladder [O. This field theory model has essentially 



the same structure as ( |4.5|) . However, in order to discuss the hole-doped case, the use 
of CP 1 variables is indispensable. 

The nonlinear-cr model with the #-term has been studied both analytically and 
numerically as an effective field theory of the AF Heisenberg models[jl3|. Properties of 
the ground state and excitations strongly depend on the value of 8, i.e., the coefficient 
of the topological term, / drl z (6) = iOQ. It is expected that the model exhibits a 
phase transition at 9 — (2m + with m = integer as 6 varies. Here we explain this 
transition from the view point of gauge theory. It is well known that a composite 
gauge field A^ appears in the CP 7 ^ 1 model, A^ = ilzd^z), which transforms as 
Ap — > Ap — d^cf) under z — ► e^z. For the case of 9 = (mod 2tt), the z-boson 
becomes massive due to the CP 1 constraint 113] and, after the ^-integration, the 



composite gauge boson acquires the Maxwell term and so becomes dynamical. In the 
(1 + 1) dimensions, the Maxwell gauge theory has only one phase, i.e., the confinement 
phase with a linear confining potential V(r) oc r between a pair of charged particles. 
Therefore possible excitations in the CP 1 model are the gauge-invariant triplet boson 
n = (zaz). These results are obtained for the CP^ -1 model by the 1/iV-expansion 
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I5| , but expected to be correct for the N = 2 case also. 

For 9 = 7i (mod 2tt), it is expected that the behavior of the dynamical gauge 
field is quite different from the 9 = case. In the present CP 1 model, the #-term is 
rewritten in terms of the gauge field (in the real-time formalism) as 

9 



f dtl z {9) =i— [ dxdtE, (4.6) 

J Z7T J 



where E is the electric field. 



Coleman JT6[ gave the following semi-classical argument on the gauge dynamics 
of QED 2 with the above 0-term. As E = d x A Q , it is obvious that the effect of the 
#-term is interpreted as putting ±9/(2tt) charges at the spatial infinities x = ±oo. In 
the case of 9 = 0, an electric flux appears through the Gauss' law between a pair of 
oppositely charged sources (say, an electron and a posotron), giving rise to a linear 
potential which confines an electron and a positorn in one-spatial dimension. When 
a pair of ±9/2tt charges are put in the spatial infinities, this confinement picture is 
not changed till the magnitude of charges increases up to 9/2n = 1/2. In the case 
9/2tt = 1/2, the "ground state" supports an electric flux of magnitude 1/2 lying along 
the entire space. As a pair of electron and positron are put into the system at x±,X2, 
there appear step-function-like jumps in the electric flux at Xi,x 2 - The electric field 
is E(x) = 1/2 for x < x\ and x 2 < x and E{x) = 1/2 — 1 = —1/2, but its magnitude 
is still 1/2 for all points[|Tj]]. Therefore, the energy, proportional to E(x) 2 , does not 



change as x±,x 2 are varied, so there is no confining force between charges. Thus, 
at 9/2tt = 1/2, the system exhibits a phase transition from the confinement phase 
to the deconfinement phase. Similar behavior is expected also for the CP^ -1 model 
with the #-term. The unsolved problem is whether the phase transiton at 9 = ir 
(mod 27r) is of first order or of second order. It may depend on the magnitude of the 



coupling constant [13]. However, it is known [|l8j that the S = 1/2 AF Heisenberg 
chain has gapless modes and spin-spin correlation functions have a power-law decay. 
These low-energy properties are described by the k = 1 Wess-Zumino-Witten model. 
Therefore it is correct that the CP 1 model with a #-term that corresponds to the 

11 



S = 1/2 AF Heisenberg chain must have a second-order phase transition at 9 = tt. 

The above discussion on the composite gauge boson reveals the essential difference 
between the spin chain and the spin ladder. As we showed, one has 9 = rr for the 
chain, which leads to the deconfinement phase of composite gauge boson, whereas one 
has 9 = 2ir for the two-leg ladder, which leads to the confinement phase. We expect 
that also for the hole-doped case the #-term is ineffective in the ladder t-J model and 
one may ignore its existence. Detailed study of the dynamics of the composite gauge 
boson and quasi-excitations in the ladder t-J model with doped holes will be given in 
the following section. 

We have considered the case J > J' so far. Similar results are obtained for the 
case J' > J. The continuum limit of the spin part is again the CP 1 model with 9 = 
in (|4.5|) , but the effective coupling constant g 2 is given as [19| 



92 = 1/5 + Tr < 4 ' 7 > 



and v z = y2J(J + y)a. It is interesting to notice that for 9 = (mod 2n) the 
QpJV-i model is asymptotically free and g 2 becomes large at low energies. This fact 
and Eqs.( |4.5| ) and ( |4.7| ) suggest J' / J — > +00 for the low-energy limit, as it is expected 



from the appearance of the energy gap for any finite value of J' (7^ 0) in the ladder 
model. (The two points J' = and J' = 00 may be fixed points of renormalization 
group. ) 

The effect of doped holes on the dynamics of spins is examined by integrating out 
the hole field in the effective model. The relevant terms come from K + T\ in ( |3.16|) , 



exp 



J drAA z ] = J [Dip] [Drjl exp [ / dr(K + T\)] . (4.8) 



To this end, the hole-hopping expansion is quite useful and reliable especially at 
sufficiently high T [XT]. In the hopping expansion, the bare propagator of holes is 



obtained from the .fT-term (|3.16|) as follows; 
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e 



-m(n— ti) 



9{n ~ r 2 ) - e-P m 6(r 2 - n 



1 + e-P m 
G v (r 1 -r 2 ) = (ViA^ViA^)) 

= -G;(r 2 -n), (4.9) 

where we employ the temporal gauge A = for simplicity. The hole concentration 
is expressed as 

<^> - 5 = YT ^- (4-10) 
From ( |4.8| ) and ( |4.9| ), we obtain AA 2 for high T as follows; 

f drAA z = j drdr'J2tX{r)c v {T') ■ (^ (t)t? Vv(t')> 

~ 5(1 - 5)/? / drX;^(r)c v (r). (4.11) 

It is straightforward to obtain the continuum limit of AA Z (4.11) as 



AA Z = J dx\C T a- x D T zD r z + C x aD x zD x z] , (4.12) 

where 

Thus the effect of hole hoppings is incorporated in the form of renormalization of the 
relativistic CP 1 model but the #-term is not generated. This renormalization increases 
the effective coupling g 2 , hence increases the spin gap. This is expected since hoppings 
of holes should reduce the spin-spin correlations. However, as T is lowered, AA Z in 
Q4.12j ) becomes to overestimate the hole-hopping effects, not only because one needs 
to include higher-order effect of hopping expansion, but because of correlations among 
holes. As we explained before, there exists attractive force between holes sitting on 
the NN sites, as T term in ( |3.16| ) indicates. This effective interaction between holes 
generates correlations among NN holes to hinder the single-hole hoppings. In some 
region, it may give rise to hole pairings and superconductivity as discussed in Sect. 6. 



13 



5 Dynamics of composite gauge bosons and quasi- 
excitations in the confinment phase 



In this section we shall study the gauge dynamics of the t-J model on a ladder. In 
Ref. we studied the CSS of the t-J model on a 2D lattice both in the slave-boson 
and SF representations, and calculated the transition temperature T css (<5) of the 
confinement-deconfmement (CD) transition. It was shown that the CSS takes place 
at low T's below Tcss, being compatible with experiments. The method presented in 
Sect. 4 of Ref. J7| can be applied in a straightforward manner to the study of the CSS 
in the present ladder system. Therefore we present below the main steps and results. 
The reader who wants to know more details should refer Ref. . 

The effective lattice model of Sect. 4 treats even and odd sites asymmetrically by 
integrating out the odd-site spins. To make the calculations and presentation below 
simpler and more transparent, we introduce another effective lattice model, which 
treats even and odd sites in a symmetric manner. The new symmetric model is very 
naturally obtained from the asymmetric model by adding odd-site spin variables so 
that it recovers the even-odd lattice symmetry. Both the previous asymmetric model 
and the present symmetric model have the same naive continuum limit in the spin 
part, i.e., the CP 1 non-linear sigma model. Since the relation between the asymmetric 
model and the symmetric model are so intimate, these two models must fall into the 
same universality class; in particular the qualitative result on the CSS derived below 
should apply also to the asymmetric model. 

To avoid confusion with the effective model in Sect. 3 and to make the expressions 
more transparent, we change some notations; For example we use x ( and y) to denote 
sites of the ladder. The partition function of the symmetric model has the following 
path-integral representation; 

Z = J [dz] [d(] exp{J Q drA) , 
A = A z + A ( , 
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A z — j,^|D r ^ z | +^ _ \ z y z x\ > 

4J * (xy) A 

A c = -J2( X (D T + m x )( x 

X 

txyipxyCyCx C-xyCxCy) i (^-l) 

{xy) 

where J2(xy) denotes summation over NN pairs (xy) on the ladder [p0[| . We write 
J = (2J+ J')/4, and (t^) implies J(t) for vertical pairs and J'(t') for horizontal 
pairs. (Grassmann number) denotes 

f ^ x hole at even site, 

I 77^ anti-hole at odd site. 
Due to this rewriting, the fermion mass becomes staggered, i.e., m x = m for even 
sites and m x = —m for odd sites. 

The spin part has been simplified to the CP 1 lattice model by keeping in A\ only 
the first term. The hole part is simplified by (i) ignoring T and T 2 focusing on the 
region out of the superconducting state, and (2) modifying b u and c v as 

bxy Z x Zy -\- -—^(ZyZ x ), 

Cxy Z x Zy -\- -—z(z x Zy) . (5-3) 

The action is invariant under the local U(l) gauge transformation 

z xa ► exp^i^j.) z xa , 
( x -> exp(i9 x )( x , (5.4) 

with time-independent rotation angles 9 X . 

Next, we introduce the following four composite gauge variables on the link (x, y); 



B x y 


{Z X Zy) 


-e* 


?* f> Jdy 

J-'xy^- i 


D X y — 


(Z X Zy) 


— > e~ 


-iO x /-) JOy 
±y X y^ 


F = 

J- xy — 


CxCy 


-e* 


1x p pidy 
- 1 xy^ i 


G xy = 


CxCy 


— > e _ 


-iOxr< pidy 

KJ X yK^ 



(5.5) 
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where their transformation laws are also indicated. We enforced these relations 
strictly via delta functions like S(B xy — z x z y ), etc., which amounts to insert into 
Z the following identity; 



J] J [dP] [dQ] [dR] [dS] J [dB\ [dD] [dF] [dG] 

{xy)r 

exp J dr i[P(B - (z x z y )) + Q(D - (z x z y )) 



(xy) 

+ R{F-Uy) + S{G-U y ) + h.c.} ixy) . (5.6) 
Then by taking the temporal gauge, we have 

Z = II J[dP\[dQ\[dR\[dS\ J[dB\[dD\[dF\[dG\ 

(xy)r 

x I C (R, S)I Z (P, Q, F) exp J dr ]T A xy , 
I C (R, S) = J K] exp J dr[- ]T C*(Ar + m x )( x 



i y ' (P-xyCxCy ~l~ ^xyCxCy h.C.)], 

(xy) 



I Z (P, Q,F) = f [dz] exp f dr[-i ^ l**| 2 + ~T Y^xyF xy z y z x - h.c.) 

J J 4-/ T- A, J f \ 



4J x 2J ( xy ) 

2 ^ ^ ( PxyZx^y ~t~ Q xy Z x Z y ~\~ /i.C.)], 

(xy) 



A xy = ^i\D xy \ 2 + t xy (FB + h.c.) 



2 



xy 



+ i{PB + QD + RF + SG + h.c.) xy (5.7) 

Next, we integrate over ( x and z x by the hopping expansion; an expansion in 
powers of gauge fields. To do this we need the following on-site Green functions; 

((x(Tl)( y (T 2 )} = 5 xy G x (T 1 -T 2 ), 

G x (r) = \ ± ^nr) 



P n^oo iuJ n + m x 

sgn(m x ) E(-) L exp[-m x (r + PL)]9[m x (r + (3L)] 



L 

exp(-m x T) 
1 + exp(-m x (3) 
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[9(r)-e(-r)eM-m x P)i (5.8) 



for \t\ < (3, and 

(z^ir^Zysfa)) = 6 xy 5 aa >G(T 1 -T 2 ) 



P n=-oo UJ n + (X 

2 T 

= ^exp[- v ^(r + /3^)] 

A' 

^<5(r) for » 1 

CT V (5.9) 

g for « 1 

Value of a can be estimated in the large- N approximation, etc. ||, but it is not 

necessary to derive the main conclusion below as long as it is finite. 

Up to the second order in the gauge variables, we obtain 



[dB] [dD] [dF] [dG] [dP] [dQ] [dR] [dS] exp / dr £ [A (z + 



+ ^L\D\ 2 + t xy (FB + h.c 



2 

+ i{PB + QD + RF + SG + H.c.)] xv , 
A (z = - Cl p\Ro\ 2 ~ c 2 \S\ 2 + c 3 \S\ 2 - c A \F\ 2 - c 5 \F\ 2 - c 6 \P\ 2 

+ c 7 \P\ 2 - c s \Q\ 2 + c 9 \Q\ 2 - ic w {FP - h.c.) (5.10) 
where the coefficients ci 2 ,...,io ar e given by 



ci = g exp(-/3|m|), g = (1 + exp(-/3|m|) 

C2 — 77i f, C 3 — — : — , C 4 — — 7-p-, C 5 



-2 



2\mV 6 8|m| 3 ' 4 a 1 ^' ~° 2( x 3 / 2 ' 

J ^itJ t n^y J 

a 3/2i °7 = Cg = -^rj^, Cio = 

1 r/3. 



C 6 - C 8 - -Jrt> C 7 - C 9 - -— , Cio - ~~~T7^ (5.11) 



i?o is the n = mode of the Fourier coefficients, R n = (3 Jq dTR(r)exp(—iuj n T). 

The integrations over P, Q, S are Gaussian, and the i?(r)-integration is done in 
Fourier components R n ; 

f n^nexpt-Ci/^ol 2 + l(3j2( R n F -n + h.C.)] 

cx exp(--|F | 2 )n^n) (5.12) 
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We further integrate over the pair fields F and G of fermions. To discuss the 
phase dynamics of residual gauge variables B and D, we set their amplitudes to 
be the constants PB, q {= \B xy \) and PD,q(= \D xy \), which can be determined by a 
straightforward mean field theory ignoring phase fluctuations of mean fields. The 
suffix q =(u,s) distinguishes whether the link (x,y) is vertical, (x(i,a),y(i + 1, o)) or 
horizontal (x(i, l),y(i, 2))p0|. Thus we write 

B xy = PB,qUxy,U xy G U(l), 

D xy = p D ,qV xy ,V xy eU(l). (5.13) 
Due to the compliteness condition, z a z a > + z a z a ' = b aa ' , we have the relation, 

Pl,q+Pl,q=^ (5-14) 

After some calculations, we reach the following effective lattice gauge theory; 
Z = JidU}[dV}exp(JdTY,A uv ) , 

AuV — —CD,q\U xy \ 2 — CB,q\V xy \ 2 
°D,q = \-a)PD& °B,q = [~a + /o_N2 1 , o[ „ \Pb,o: (5-15) 



Following the Polyakov-Susskind method |p2TJ for studying CD transitions in lattice 
gauge theories, this system can be mapped to the anisotropic XY model [^]; 

Z X y = J II ex PE J hq cos(a y - a x ) + J 2,q cos(a y + a x )], 
x (xy) (xy) 

Jl,q = P~ l CD,q, h,q = (3~ X C B ,q- (5.16) 

The J2 term, coming from \V\ 2 term, expresses an anisotropy of XY spin couplings, 
and reduces the global spin symmetry from U(l) (for J2 = 0) down to Z(2) (for 
J2 7^ 0). A possible phase transition of this spin system describes the CD phase 
transition of the original t-J model. For two and higher dimensional lattices, this 
spin system is known to have an order-disorder transition at J\ + J 2 — 1 0, which 
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leads to existence of a finite transition temperature Tcss- For T > Tcss, the XY 
spin correlations decay exponentially, /(|x|) = (S x Sq) ~ exp(— a\x\), which implies 
the potential energy among two gauge charges, W(|x|) = — ln/(|x|) (see Ref.f21~H for 
this relation), to be a linear-rising confining potential, iy(|x|) ~ a\x\. So the gauge 
dynamics realizes here in the confining phase. For T < Tcss, there is a long-range 
order, = const + exp(— c\x\), which implies a short-range potential ~ 

exp(— c\x\). So the gauge dynamics realizes in the deconfining phase here||22||. 

In the present ladder system, the XY model Q5.16| ) is put on a coupled two one- 
dimensional chains and its universality class is classified into that of a single one- 
dimensional chain, the partition function of which is given by 

Z c xv m = /lly !ex P[ J iE cos (« !+ i-« i ) + -h E cos («m ( 5 - 17 ) 

i i i 

One can solve its spin correlations exactly, which decay exponentially, = 
A exp(— a\x\). The coefficient a is finite as long as are finite. Thus, we con- 
clude that the gauge dynamics of the t-J ladder is always realized in the confining 
phase except for T = at which J 12 diverge. This is the main conclusion of this 
section. The case of T = is discussed later on. 

Here we comment that the above conclusion is based on the hopping expansion. 
If a, the mass of z, vanishes, the hopping expansion cannot be applicable. Our 
assumption a > is supported by the large-iV analysis of the CP^V model in 
one spatial dimension 0. This is in strong contrast with the case of a single chain. 
There, the topological term appears and it generates massless excitations, preventing 
us from making use of the hopping expansion and obtaining an effective lattice gauge 
theory. However, as discussed in Sect. 4, we know that the spin chain system is 
realized in the deconfinement phase. We expect that the doped chain still supports 
this deconfinement phase, although a solid analysis is required. (The special case of 
supersymmetry J/t = 2 can be solved exactly by the Bethe anzatz and its result 
supports the CSS.) 
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How about the gauge dynamics on a ladder at T = ? Since the system ( |5.17| ) 
develops a long-range order, one may conclude that the gauge system is in the decon- 
finement phase. However, the expression (|5.17| ) itself is not adequate at T = 0. This 
can bee seen as follows. At T = 0, m = as one can see from ( |4.10| ). So the hopping 
expansion with respect to fermions is not applicable. In this case, we have a sepa- 
rate argument using the massless Schwinger model f23|, the massless QED 2 without 
topological terms, which is obviously a relevant model for present argument. Here 
it is known [53 that the physical spectrum consists of a massive boson, a composite 
of fermion and antifermion. Since this is a gauge-invariant neutral object, its gauge 
dynamics is in the confinement phase. From this fact, we conclude that the t-J model 
on a ladder realizes its gauge dynamics in the confinement phase not only for T > 
but also for T = 0. 

How is the effect of weak but finite couplings among NN ladders in a plane and 
couplings among interlayers ? They put the system into three dimensional one. The 
calculation given above can be redone to reach the same expression ( 5.16 ) for the XY 
spin model, but the sum over (xy) is extended to the entire lattice. One concludes 
that there is a finite critical temperature Tcss below which the deconfinement phase 
appears. The value Tcss is small and approaches zero as these three-dimensional J 
and t couplings tend to vanish. One can determine Tcss using the 3D XY model and 
the MF values of ps, g , etc. 

What are the quasi-excitations in the confinement phase ? In this phase, due 
to the confining force, only charge-neutral objects can be physical excitations. The 
typical examples are the original electrons z x ip x which carry the electric charge — e 
and spin S = 1/2, the charge density fluctuation ^) x ^ x , and the spin triplet excitation 
z x az x , i.e., charge-neutral magnons of S = 1. Together with these local combinations, 
one may conceive nonlocal gauge-invariant combinations like z x U xy tp y , etc. Real low- 
energy excitations are linear combinations of these states with definite electric charge 
and spin. 



20 



Finally, we note that the analysis in this section can be repeated when one includes 
the To and T 2 terms into the symmetric model. These terms are necessary to generate 
the superconducting phase as we shall see in the following section. The corresponding 
calculations of Tqss f° r the 2D t-J model are found in Ref. [0]. After tracing these 
calculations, one reaches the same effective XY spin model of (|5.16 ) with the modified 
coefficients Ji j2 . Thus we obtain the conclusion that the confinement phase obtained 
above continues to exist even when the system is superconducting. 



6 Superconductivity 

As we saw in Sect. 4, there are effective interactions between holes which appear as a 
result of SRAFO. In terms of the composite hole-pair field M ou , T and T 2 terms in 
( |3.16j ) are rewritten as 

T 2 = - 2 T ^a-M W M W ,. (6.1) 

It is obvious that T induces a condensation of the hole-pairs, 

(M ou ) ^ 0, (6.2) 

and T 2 favors the "d-wave" symmerty, 

"£t v (M w ) = 0. (6.3) 

Fort' = t, Eq.Q means (M os ) = -2(M ou ) = -2(M od ) if we assume (M 0iU ) = (M 0jd ) 
and (M 0>s (u,d)) are independent of the position o |24| . 

It is straightforward to introduce a hole-pair field A ou by a Hubbard-Stratonovich 
transformation. We shall obtain a Ginzburg-Landau (GL) theory of the hole-pair 
field by integrating out the hole field. To this end we take the temporal gauge A = 
for calculational simplicity^]. The T term is rewritten by A ou as 



exp( J drT ) = exp ( J dr ^ y (M ou A ov + M ov A ov - A ov A ov )) , (6.4) 
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up to an irrelevant constant. From this, one gets the relation, 

(M ou ) = A ov . (6.5) 

The GL potential energy, Vgl(A) of A ou (that appears in the action integral / drA 
in the form of — /3Vql ), is obtained by the hole-hopping expansion as 



V Gh {A) = V + V 1 + V 2 + O(A 4 ), 

— ~ JuA ou A ou , 

* ov 

VX = ~\{b 2 f dTdT'Y,Jl&ov&ov-{M ov {T)M ov {T')), 

1 1 2 1 r 
V 2 = p (-) j - J//2 J drdr'dr" g J v J v ,t v t v ,A ov A ov 

x{M 0V {t)M 0V ,{t)M 0V {t')M 0V :{t'')), (6.6) 

where A ov are assumed to be time- independent. The above correlation functions of 
M 01/ 's are calculated by the hole propagator (|4.9Q and the expectation value of z-pair, 

(zz ou ) = D u . (6.7) 



The actual value of D v should be obtained by the MFT|26[, but here we only needs 



that it is nonvanishing due to SRAFO. After some calculation, we obtain 
8m ^ 

v 2 = - j^j, J m2 g (6 - 8) 

We are interested in the case of low hole concentrations i5 C 1 where we have the 
expression, 

m ~ — - — , (6.9) 

P 



which comes from ( 4.10| ). Since m is large for 8 <C 1, we have neglected the terms of 
0(m~ 3 ) in the above. 
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Let us find a solution that minimizes Vgl(A). It is easily found by minimizing 
Vo + V\ and V2 separately. By writing A ot , = e^A , we have 

Vo + V, = N Vmk xC 2 \A\ 2 

C2 = Y,{-\j, + ^\DA 2 Jl)el. (6.10) 

V2 is minimized by configurations satisfying 

Y,t„JvD v A v = 0. (6.11) 



Eq.( |6.10D implies that the order parameter A devevlops for C2 < 0. So the critical 



temperature T c is calculated from C2 = as 

T = Wei 1 f , 

T^Juel 4|ln5| 4|ln5|- 1 ' J 

The last expression is for the case of J' = J , t' = t and ^-independent e v and D v = D. 
In this case, Eq. (|6.11| ) gives rise to the <i-wave solution [ pjj] , 



A os = -2A ou = -2A od . (6.13) 

The MFT using a Gutzwiller renormalization of matrix elements || predicts similar 
behavior of the superconductivity correlations. If we take into account the fluctuations 
of the order parameter of the superconductivity A ou , weak but finite three-dimensional 
inter-ladder interactions are required to realize a genuine superconducting phase. 
From our discussion given so far, it is obvious that the existence of the SRAFO is 
essential for appearance of the superconductivity and its d-wave type symmetry. 



7 Discussion 

In this paper, we studied the t-J model on the ladder by using the SF-CP 1 formalism. 
We obtained the effective low-energy model by integrating out the half of the CP 1 
variables. At low T's, there exists the SRAFO which makes it easy to integrate 
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over the CP 1 variables. The low-energy effective model shows why the low-energy 
excitations are different in the ladder and the chain systems, reflecting the coefficient 
of the topological term in the CP 1 model. In the ladder system, the confinement 
mechanism works at any T and therefore excitations are spin triplet and spinon- 
holon bound state. On the other hand, in the chain case, the full CSS takes place 
because of the ^-term of the composite gauge boson. Similar analysis works for the 
2D t-J model, which is studied in the previous papers. We stressed that the gauge- 
theoretical point of view, especially the CD transition, is quite useful and universally 
applicable for such problems of separation of degrees of freedom. 

For even-number-leg ladder Heisenberg models, it is expected that an energy gap 
appears for spin excitations for any finite value of J'. The continuum field theory 
model is the CP 1 model which is known to be asymptotically free. Therefore the 
coupling constant becomes large at low energies. Our derivation shows that this 
means J'/ J — > +00 for low-energy limit as expected. 

The effective model also shows that there are effective attractive interactions be- 
tween holes in the SRAFO background. By these interactions, the superconducting 
phase is possible under weak but finite 3D inter-ladder interactions. 

In conclusion, we observed that the SF-CP 1 formalism and the gauge-theortical 
method provide a natural and coherent way to understand various important prop- 
erties of the t-J model on a ladder. 
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Figure Caption 



Fig.l Layout of a ladder. Two vertical lines denotes two legs, a = 1 and a = 2. 
A site on the ladder is labelled by a set of suffices (i,a), where % counts the site 
along each leg. The filled circles denote odd sites, and the open circles denote even 
sites. For a fixed odd site o, there are three nearest neighbor even sites, (o, v) with 
v = u (up), v — d (down), v — s (side). For example, for o = (i, 1), (o,u) = 
(i + 1,1), (o, d) = (i — 1,1), (o, s) = (i,2) as illustrated in the Figure. For o = 
(j,2), (o,u) = (j + 1,2), (o,d) = (j - 1,2), (o,s) = (j, 1). The hopping amplitude 
and the exchange coupling are t, J, respectively, for vertical links (i, a), (i + 1, a), and 
t', J' for horizontal links (i, 1), (i, 2). 
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